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Abstract. By definition a separable state has the form 5Zf =1 Ai® 
Bi, where < Ai,Bi for each i. In this paper we consider the class 
of states which admit such a decomposition with B\ , . . . , B p having 
independent images. We give a simple intrinsic characterization of 
this class of states, and starting with a density matrix in this class, 
describe a procedure to find such a decomposition with B\ , . . . , B p 
having independent images, and A±,...,A p being distinct with 
unit trace. Such a decomposition is unique, and we relate this to 
the facial structure of the set of separable states. 

A special subclass of such separable states are those for which 
the rank of the matrix matches one marginal rank. Such states have 
arisen in previous studies of separability (e.g., they are known to be 
a class for which the PPT condition is equivalent to separability). 

The states investigated also include a class that corresponds (un- 
der the Choi-Jamiolkowski isomorphism) to the quantum channels 
called quantum-classical and classical-quantum by Holevo. 



Introduction 

A state on M m g) M n is separable if it is a convex combination of 
product states. States that are not separable are said to be entangled 
and are of substantial interest in quantum information theory since 
entanglement is at the heart of many applications. Some useful nec- 
essary conditions are known for separability, e.g., the PPT condition, 
by which a separable state must have positive partial transpose [Hj. 
There also are some necessary and sufficient conditions, e.g. [9], which 
however are difficult to apply. Thus it would be of great interest to 
find a practical test for separability, at least for a significant class of 
states. 

Closely related to this is the goal of finding a procedure to decom- 
pose interesting classes of separable states into a convex combination 
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of product states. Such a procedure would not only shed light on sep- 
arable states, but would provide a separability test for that class. 

We will identify states with their associated density matrix, and also 
consider unnormalized states, which are then associated with positive 
semi-definite matrices. (We will abbreviate "positive semi-definite" to 
simply "positive" hereafter.) Thus a density matrix T e M m <g> M n is 
separable if it admits a representation 



where each A y and _B 7 is positive. Such a density matrix T represents 
a mixed state on a bipartite quantum system composed of two sub- 
systems, the A-system and the 5-system, associated with M m and M n 
respectively. 

In our previous paper [2], the authors studied separable states with 
such a representation with each and B 1 rank one, with the require- 
ment that Bi, . . . , B p be projections onto linearly independent vectors. 
This class of states turns out to be same as the set of separable states T 
with the property that T and the marginal state tr^ T = Tg obtained 
by tracing out over the A-system have the same rank, cf. Lemma [15j 
The equivalence of these two formulations was established for states on 
M2 ® M n in [13], and then in complete generality in [101 Lemma 6, and 
proof of Thm. 1], where it was also shown that for states satisfying 
this rank requirement, the PPT condition is equivalent to separability. 
(An alternate proof of the equivalence of these rank and independence 
conditions was given in [THl Lemma 13].) In [TU] the authors also gave a 
procedure for decomposing such states into a convex sum of pure prod- 
uct states, based on an inductive argument for finding a certain kind of 
product basis, and then a reduction to a block matrix whose blocks are 
normal and commute. In this paper we also make use of a reduction 
to this type of matrix. The existence of special families of commut- 
ing normal matrices played an important role in the investigation of 
separability in [17J as well. 

The current paper investigates separable states for which there is no 
rank restriction, but admitting a representation (pQ) in which Bi, . . . , B p 
have independent images. We call such states S-independent, and give 
an intrinsic way to determine if a state falls in this category (Theorem 
[6]). We show that without knowing an explicit decomposition to begin 
with, there is a canonical way to locally filter T to yield a state T which 
admits a representation (pQ) in which B\, . . . ,B P are orthogonal. (Of 
course, there is nothing special about the 5-system compared to the 
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A-system, and all results in this paper are valid with the roles of the 
A and B systems interchanged.) 

This is then used to give a canonical form for T, and to find a decom- 
position of T of the form (HI), cf. Theorem [6j This decomposition can 
be chosen so that Ai, . . . , A p are distinct and have unit trace, and in 
that case the representation is unique. It is then simple to decompose 
further to get a representation of T as a convex combination of pure 
product states (i.e., of density matrices where each is the projection 
onto the span of a product vector), and we describe when this decom- 
position is unique (Theorem UM) Finally, we show in Theorem HT1 that 
if a state has a representation (JTJ) with the images of the A 1 disjoint 
and the images of the B 1 independent, then the face of the space S 
of separable states that is generated by this state is the direct convex 
sum of separable state spaces of lower dimension. 

The density matrices investigated here are closely related to inter- 
esting classes of completely positive maps. A completely positive map 
$ : M m — > M n is entanglement breaking if (/ (g) $)(r) is separable for 
all positive T, cf. [THEE]. The Choi-Jamiolkowski isomorphism [3] IT2] 
is a linear isomorphism under which completely positive maps corre- 
spond to positive matrices. Under this correspondence, entanglement 
breaking maps correspond to separable matrices, so the results of this 
paper on convex decompositions of a class of separable states then can 
be transferred to give information about decompositions and identifi- 
cation of the corresponding class of entanglement breaking maps. 

In particular, there are two important classes of entanglement break- 
ing maps (quantum-classical channels and classical-quantum channels) 
that have Choi matrices in the class of separable states investigated in 
the current paper. These classes were originally singled out by Holevo 
[8], and further investigated as special cases of entanglement breaking 
maps by Horodecki, Shor, and Ruskai in [TTj [13] . These are shown 
in Theorem [8] to be special cases of the classes of A-orthogonal and 
5-orthogonal density matrices, which play a key role in the current 
paper. Theorem H] and Theorem [8] together provide an intrinsic way 
to identify such quantum channels without knowing an explicit Kraus 
decomposition ahead of time, as well as giving a procedure to find a 
Kraus decomposition of the appropriate form. 
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A CLASS OF SEPARABLE DENSITY MATRICES 

Definition. Subspaces V\, . . . , V p of a vector space are independent if 
their sum is a direct sum. This is equivalent to the implication 

p 

x 7 = with x 7 G V y for 1 < 7 < p ==>- all x 7 = 0. 

7=1 

We now define the central class of separable density matrices that 
we will investigate. Later in Theorem [5J we will give an intrinsic char- 
acterization of this class. 

Definition. A density matrix T G M m <g> M n is B-independent if T 
admits a decomposition 

p 

(2) T = A i ® B i 

7=1 

where < Ay, 5 7 for 1 < 7 < p, with the images of Bi, . . . , B p inde- 
pendent. 

Example. Let xi, . . . ,x p G C' m and yi, . . . ,y p E C n be unit vectors, with 
yi, . . . ,y p linearly independent, and < Ai, . . . , X p with A 7 = 1. Let 
T be the convex combination 

p 

(3) T = A 7 ^ 7 ® 

7=1 

where for a unit vector z, P 2 denotes the projection onto Cz. Then T is 
5- independent. The uniqueness of such decompositions, and the facial 
structure of faces of the separable state space generated by such states, 
were investigated by the current authors in [2]. As was discussed in 
the introduction, such states played an important role in [10] and also 
appeared in [15] . 

We will return to the subject of £>-independent states after develop- 
ing some necessary results. 

One sided local filtering 

Definition. A linear map $ : — > M,i is a filter if there is a positive 
A G M d such that $(X) = AX A. (We do not require that A be 
invertible.) A map $ : M m Cg) M n — > M m (g) M„ is a /oca/ /i/ter if there 
are positive A, B such that $(X) = (A <S> B)X(A <g> B). 

Applications of filtering, e.g., to distillation of entanglement, date 
back at least to [3j HJ Ej. It is well known that we can apply a local 
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filter to any density matrix to arrange for one or the other partial trace 
to be a projection, as we now describe. 

Definition. If A > 0, then A* denotes the Penrose pseudo-inverse of 
A, i.e., the unique positive matrix which is zero on (im A) 1 - and satisfies 
A* A = AA# = Pa, where Pa is the projection onto the image of A. 
If a spectral decomposition of A is A = £\ AjPj with all A, > 0, then 

A * = J2i K lp i- 

Definition. We write tr# and tr^ for the partial trace maps on M m £g> 
M n , and if T e M m £g> M n then we write Tb = tr^ T and Ta — tr# T. 

Definition. Let < T e M m Cg> M n . Then we denote by T the matrix 

(4) T = (I <g> ((T B ) # ) 1 / 2 )T(J ® ((Tb)*) 1 / 2 ). 

We view the pair (T,T B ) as partitioning information about T into 
a state Tg that contains information about T on the subsystem B, 
and another part T which contains information about T relating to the 
system A as well as the interaction between A and B systems. 

We will show later in this section that T can be recovered from the 
pair (T, Tg). First, we discuss various facts about partial traces and 
filters which we need. 

If T = £V <S> G M m eg) M n , then by definition 

(5) tr A T = J2 Ta and tr B T = J2 ir ( T v) E H- 

i ij 

It is well-known that the partial trace maps are positive maps. We 
now show that they are also faithful, i.e., if T > and either partial 
trace of T is zero, then T is zero. (We expect the following is well- 
known, but we have included it here for lack of an explicit reference.) 

Lemma 1. The partial trace maps are faithful. 

Proof. Let < T = £V <g> T {j e M m ®M n . We first show 

(6) T = <^ Ta = for 1 < % < m. 

If T = 0, then clearly all Tj = 0, so in particular all Ta — 0. Conversely, 
suppose all T, are zero. Then (En <g> I)T(Eu ® I) — En <g> Ta — for 
all i. Then for 1 < i < m, 

((E ii ®I)T 1 ' 2 )((E ii ®I)T 1 ' 2 )* = 

so <g> J)T 1 / 2 = 0. Adding gives 

= J2( E H ® J ) Tl/2 = ( J ® J ) Tl/2 = 
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and hence T = 0. 

Now we are ready to prove faithfulness of the partial traces. By (jSJ), 
if tre(T) = 0, then in particular tr(T^) = for each i. Since < T, 
then < Tjj for each i, so tr(T^) = implies = for all % and thus 
T = 0. 

On the other hand, tr^(T) = implies J^i^a = 0, and by positivity 
of each Ta, we again have T„ = for each i, and thus T = 0. □ 

We next review some useful facts about projections and images. (For 
additional background, cf. [H Chapter 3].) If A = A* and P is a 
projection, then 

(7) imicimP PAP = A. 

(Indeed, if imA C imP, then PA = A, so taking adjoints and using 
A* = A gives A = AP. Then PAP = P(AP) = PA = A. The 
converse implication is clear.) 

If £ 6 M r is a projection, then we write E' — I — E, where / is 
the identity in M r . Note that for E a projection in M m , (E <g> I)' = 
(J ® /) — {E ® I) = E' £g> I. For any projection R and positive operator 
T we have 

(8) RTR = T R'TR' = 0, 

cf., e.g., [U Lemma 2.20]. 

Finally, we observe that if < Ax, A%, . . . , A p , then 

(9) imJ2 A i = J2 imA i- 

i i 

Indeed, for 1 < j < p we have Aj < Yli A i so ker^VA; C kerAj. 
Taking orthogonal complements shows imA, C im^V A ; which im- 
plies ^2 j im Aj C im^\A;. The opposite containment is evident, so 
follows. 

The following result is clear for separable T, but requires a little 
more work for general T. 

Lemma 2. If < T 6 M m ® M n? t/te minimal product subspace con- 
taining the image of T is im Ta <8> im • In particular, if Pb is the 
projection onto the image of Tb, then (I cg> Pb)T(I £g> Pb) = T . 

Proof. Let V C C m and W C C n be subspaces, and let the corre- 
sponding projections be P and Q- Then by ([7]), imT C V <E> W iff 
{P®Q)T(P®Q) = T. 

Note (P g) Q)T(P ® Q) = T is equivalent to the combination of 
(P <g> I)T(P ® I) =T and (/ <g> Q)T(7 ®Q)=T. Thus it suffices to 
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show that 

(10) (P<g> I)T(P®I) = T imP D im(T A ) 

together with the corresponding statement for Tb- Since the proof of 
the statements for T4 and Tb are essentially the same, we just will 
prove the statement for T4. 

We will make use of the following identity valid for all T e M m <g> M n 
and all X e M n : 

(11) \,t b {X®I)T(X®I)=X{\z b T)X. 
Thus 

(P ® I)T{P ®I) = T ^ (P' <g> /)(T)(P' ® /) = by © 

tr B ((P' (8) P)T(P' ® J)) = by Lemma □ 
^ P'(tr B T)P' = 0by (P) 
P tre TP = ti B T by © 



(12) 



imP D mitre T = im(T^) by ([7]). 



This completes the proof of (fTUj) . and hence finishes the proof of the 
lemma. □ 

The next result relates properties of T and T, and shows that T can 
be recovered from the pair (T, Tb). 

Lemma 3. Let < T e M m £g> M n and define T as in (j4j). T/ien 

(13) T = (/ ®T^ /2 )f(I ®T^ 2 ). 

T will be separable iffT is separable, and tr^T = P B (where Pb is the 
projection onto the image ofT B ). 

Proof. From the definition of T, separability of T implies that of T. 
For F = ((Tb)*) 1 / 2 we have 

tr A f = tr A (I ® F)T{I ® F) = P(tr A T)P = FT B T = P B . 

Furthermore, 

(/ <g> T l J 2 )f{I ® T l J 2 ) 

= (I <8> Tj /2 )(J ® ((T B )*) 1/2 )T(I ® T^ /2 )(J <g> ((T^)*) 1 / 2 ) 

(14) = (/ <g> P B )T(J <g> P s ) 

This would prove fflBl if we knew the range of T were contained in 
C m (g) im Pg . This follows from Lemma [2J Finally (fT3"|) shows that 
separability of T implies separability of T. □ 
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.B-ORTHOGONAL DENSITY MATRICES 

In this section we describe a canonical form for a class of positive 
matrices which we call .B-orthogonal, and which is a subclass of the 
^-independent matrices. In the following section we will apply these 
results to achieve a canonical representation for the full class of B- 
independent matrices. 

Definition. Positive matrices in M r are orthogonal if their images 
are orthogonal. A density matrix T is B-orthogonal if it admits a 
representation 

p 

(15) T = ^A y <8>£ 7 

7=1 

with < Ay, B 1 and with the {-B 7 } matrices orthogonal. Similarly we 
say T is A-orthogonal if it admits a representation ( 1T5|) with the {Ay} 
matrices orthogonal. 

Definition. {Eij} are the standard matrix units of M m . For any ma- 
trix S G M m g) M n , we denote by {Sy | 1 < i,j < m} the unique 
matrices in M n such that 

S = J2 E V® S V- 

The following gives a canonical form for .B-orthogonal matrices, and 
a readily tested necessary and sufficient condition for 5-orthogonality. 

Theorem 4. Let < T = J2ij E ij ® e M m ® M n . Then the 
following are equivalent. 

(i) T zs B-orthogonal. 
(iij ^4.// 

are normal and mutually commute. 

Furthermore, if T is B-orthogonal, then T admits a unique represen- 
tation 

p 

(16) T = ^A y ®Q 7 , 

7=1 

with Qi, . . . , Q p orthogonal projections, and Ai, . . . ,A P distinct nonzero 
positive matrices. 

The projections Qi, ■ ■ ■ ,Q P will be the projections onto the joint eigenspaces 
of {Tij} (excluding the joint zero eigenspace), and will have sum Pb 
(the projection onto the image ofTg). The matrices A^ are given by 

(17) A 7 = ti A (I®Q 1 )T(I®Q 1 ). 
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Proof, (i) ==>■ (ii). If ( |T5|) holds with Bx, . . . ,B P orthogonal, then for 
each pair of indices i,j 

v 

(18) T tj = tr A ((Eji ® I)T) = £ tr(E ji A 7 ) J B 7 . 

7=1 

Since Bx, . . . ,B p are orthogonal, then Bi, . . . ,B p commute. It follows 
that the matrices {T^- | 1 < i,j < m} commute and are normal. 

(ii) ==>- (i) Conversely, suppose {T^- | 1 < i, j < m} commute 
and are normal. Define Qx, . . . , Q p to be the projections onto the joint 
eigenspaces (for non zero eigenvalues) of {T^}. For each i, j write 

v 

(19) T ij = J2^ J QT 

7=1 

Then 

t = Y! E a ® 7 o = E ^ ® ( E A: ^~ 

ij ij \7=1 



( 20 ) =E E A 7^ ®«7- 



For each 7 define A 7 = A^'E^ e M m . Then T = £ 7 ^ 7 ® <2 7 - 
For each i, j, 7 we have A 7 J = (Ay)^. Thus by the definition of the joint 
eigenspaces of {T^}, for 71 7^ 72 we must have A yi 7^ A 72 , and hence 
Ax,...,Ap are distinct. Now orthogonality of Qx ■ ■ ■ , Q p implies (TT7j) . 

Finally, we prove uniqueness. Suppose that we are given any repre- 
sentation ffl6|) of T where {<5 7 } are orthogonal projections and {Ay} 
distinct nonzero positive matrices. Then for 1 < i,j ' < m, 

7 

Then the image of each Q 1 consists of eigenvectors for for the eigen- 
values ti(EjiA^), and by distinctness of Ax, ■ ■ ■ , A p for 7! 7^ 72 there is 
some pair of indices i,j such that tx{Ej i Ar 11 ) 7^ tr(i5jjAy 2 ), so the Q 1 
are precisely the projections onto the joint eigenspaces. □ 

Remark. The condition (ii) is equivalent to the existence of an orthonor- 
mal basis of joint eigenvectors for {Tjj}, as is well known. 

A CANONICAL FORM FOR 5-INDEPENDENT MATRICES 



The following describes how to map positive matrices with inde- 
pendent images to orthogonal projections by filtering with a positive 
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matrix. We say an Hermitian matrix A e M n lives on a subspace H 
of C™ if imA C if (or equivalently, if A = on fP 1 ). 

Lemma 5. Let Xi, . . . , Xp e C n mit/i imXi, . . . , imX p independent, 
and let P be the projection on the image of^2 i X i . Then 

A=(£» # ) 1/2 

i 

is the unique positive matrix living on imP such that {AXiA | 1 < i < 
p] are orthogonal projections with sum P. 

Proof. Let A = ((T,i X i)*) 1/2 , and define Y i = AX i A for 1 < * < P- 
Then 



J2 Y > = J2 AX > A = A (j2 x ) 

i i \ i / 



A 



where P is the projection onto the image of Yli X i- By assumption, 
im Ji, . . . , imX p are independent. Since for each i, A is invertible on 
imP D imXj, and im Yj C A(imXj), then Y 1; . . . , Y p have independent 
images. Now for 1 < j < p, 



Yj — PYj — YiYj, 



and then independence of the Y's implies YiYj = for % ^ j, and 
Yf = Yj, so Yi, . . . , Y p are orthogonal projections with sum P. 

Finally, to prove uniqueness, suppose that < A , with A living on 
imP and with {AqX^Aq | 1 < i < p] projections with sum P. Then 
im/lo C imP, and 

(21) A (J2^)Ao = P, 

i 

so im Aq = im P. Multiplying (I2T!) by Af on left and right of each side 
gives J2i x i = (At?, so A = ((E,^) # ) 1/2 - ' □ 

Theorem 6. Let < T G M m ® M n . T/ie following are equivalent. 

(i) T is B -independent. 

(ii) T zs B-orthogonal. 
(hi) ^4.// 

are normal and mutually commute. 
IfT is B -independent then T admits a unique decomposition 

p 

(22) T = ^2 A-y <S> B 1 

7=1 

with < Ay, P 7 , tr Ay = 1 ; Ai, . . . , A p distinct, and B 1 , . . . ,B P inde- 
pendent. 
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Let Qi, . . . , Q p be the projections corresponding to the joint eigenspaces 
of {Tij} excluding the subspace corresponding to the zero eigenvalue. 
Then the unique decomposition (J22l) is given by 

(23) B, = (Tb) 1 / 2 <5 7 (T b ) 1 / 2 

and 



(24) A 7 = tr B (/®g 7 )T(/®g 



and the sum of the projections Q 1 will be the projection onto the image 

"/•//;. 

Proof. If T is -B-independent, then by definition there are positive ma- 
trices A\, . . . , A p and positive matrices Bi, . . . ,B P with independent 
images such that 

p 

(25) T = ^A 7 <g>-B 7 . 

7=1 

If necessary, we absorb scalar factors into the _B 7 so that tr A 1 = 1 for 
all 7, and we combine terms if necessary so that A±, . . . , A p are distinct. 
Now by the definition (jl]) of T, 

(26) f = Y j ^® ({TBffl'B^iTB)*) 1 ' 2 = A ® Q, 

7 =1 7=1 

where 

(27) Q 7 = {{Tsff^B^Ts)*) 1 / 2 . 

By Lemma O since T B = B^, then Q 1 , . . . , Q p are projections with 
sum the projection onto the image of B> 7 , and hence Q 7 = Pb- 

Thus T is 5-orthogonal. Furthermore, by the uniqueness statement 
of Theorem HI Q 7 and A 7 must be as described in that Lemma (with 

T in place of T). By (1271) . since each B 7 has range contained in the 

_____ 1/2 

range of Tg, then multiplying p7j) on both sides by T B gives (1231 . and 
ff24|) follows either from (126]) or from Theorem HI Thus we have shown 
that if T is -B-independent, then T admits a unique representation as 
specified in the theorem. 

To show that S-orthogonality of T implies ^-independence of T, we 
apply Theorem H] again. We have the representation 



T = ^ j A 1 ®Q, l 
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where A 1 and Q 1 are defined as in Theorem 0] with T in place of T. 
Note that the image of each Q 1 will be contained in the image of Tb, 
and Tb = Pb by Lemma [31 Define B\ , . . . , B p by ( 1231) . Orthogonality 
of the Q 1 implies that their images are independent. By definition, 
(Tg) 1 / 2 is invertible on the range of Tb, and imB 7 C (T B ) 1 / 2 (imQ 7 ), 
so E>i, . . . , B p have independent images and are positive. By equation 
(H3D of Lemma El 



so T is ^-independent. 

Finally, equivalence of (ii) and (iii) follows from Theorem HI 



Connections with QC and CQ quantum channels 

We will show in this section that the quantum channels known as 
classical-quantum channels and quantum-classical channels correspond 
under the Choi-Jamiolkowski isomorphism to density matrices that are 
in the classes of matrices we have called A-orthogonal or I?-orthogonal 
respectively. The remainder of this paper is independent of this section. 

Definition. Let $ : M m — > M n be a quantum channel (i.e., a com- 
pletely positive trace preserving map). If it is possible to choose < 
Fx, . . . , F q e M m , < Rx, ■ ■ ■ , R q , and tr Rk = 1 for all k such that 



fe=i 

such a representation is called a Holevo form for $. (Note that since 
$ is assumed to be trace preserving, we must have Ylk^k — 1 ■) 

The following notion is due to Holevo [8], and was further investi- 
gated in |TT] in the context of entanglement breaking maps. 

Definition. A quantum channel $ : M m — > M n is a classical- quantum 
(CQ) channel if $ admits a Holevo form fl28|) with Fx, . . . , F q rank one 
projections (necessarily with sum I m since $ is a quantum channel). 
Similarly, one says $ is a quantum-classical (QC) channel if $ admits 
a Holevo form with Rx, . . . , R q rank one projections with sum I n . 

Definition. If $ : M m — > M n is a linear map, the associated Choi 
matrix is the matrix in M m £g> M n defined by 



T = {I® T l J 2 )T{I ® T l J 2 ) = J2 A ~f® R 



□ 



(28) 
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where {i%} are the standard matrix units of M m . 

It was shown by Choi [5] that $ is completely positive iff C$ is pos- 
itive semi-definite. Note that $ will be trace preserving iff tr$(i£y) = 
Sij, or equivalently, iff tr B C$ = /. 

Lemma 7. Let F 1 , . . . , F q e M m and R 1 , . . . , R q € M n . Define $ : 
M m ^M n by 

$(X) = ^tr(F fc X)i? fe . 

k 

Then the corresponding Choi matrix is 
(29) C^ = Y, F i®Rk- 

k 

Proof. This follows from (T8J Theorem 2 and Lemma 5], or directly 
from the definition of the Choi matrix: 

c$ = E a ® $ (^y) = ^ ® tT ( E n F k)Rk 

ij ij k 

k \ ij / 
k \ ij ) 

= Y,Fl®R k , 

k 

where the final equality follows from the fact that the matrix units 
{E^} are an orthonormal basis for M m with respect to the Hilbert- 
Schmidt inner product. □ 

Theorem 8. Let < T e M m ® M n . 

(i) T is the Choi matrix for a QC channel iff T is B-orthogonal 
with tr B T = I . 

(ii) T is the Choi matrix for a CQ channel iff T is A-orthogonal 
with tiB T = I . 

Proof, (i) Let $ : M m — > M n be a QC channel with Choi matrix T. By 
definition, there is a Holevo representation (12"5|) with Ri, . . . , R n rank 
one projections with sum I n . By Lemma [7] the Choi matrix for $ is 

n 

T = J2 F l® R k- 

k=l 
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Since J2i Ri = In, then Ri, . . . , R n are orthgonal, so T is 5-orthogonal. 
Since $ is a quantum channel, then tipT = J. 

Conversely, suppose T is 5-orthogonal with tr# T = I and rank tr^ T = 
n. Since T > 0, then <3> is completely positive, and since trg T — I, 
then T is trace preserving, so T is a quantum channel. By definition 
of B-orthogonality, we can write 

p 

T = Y^A k ®B k 

k=l 

with < Ax, . . . , A p and < B\, . . . , B p with B\, . . . , B p orthogonal. 
Via its spectral decomposition, we replace each Bj by a linear com- 
bination of orthogonal rank one projections, and absorb scalar factors 
into the Aj's. Then we can write 

q 

(30) T = J2 F l®Rk 

j'=i 

with Ri,...,R q orthogonal rank one projections. Clearly q < n. If 
q < n, we can define F q+ i, . . . , F n to be zero, and choose rank one 
projections Rq+i, . . . , R n so that J2i ^» = In- Thus $ admits a Holevo 
form f)30p in which Ri, . . . , i? n are rank one projections with sum I n , 
so $ is a QC channel. 

The proof of the characterization of CQ channels is similar. 

□ 

Faces of the separable state space 

A face of a convex set C is a subset F such that if A and B are 
points in C and a convex combination tA + (1 — t)B with < t < 1 
is in F, then A and I? are in F. The intersection of faces is always a 
face, so for each point A e C there is a smallest face of (7 containing 
^4, denoted facec A 

We let K (or fQ) denote the convex set of states on Md, i.e., the 
density matrices, and S (or S mn ) denotes the convex set of separable 
states on M m ® M n . There is a canonical 1-1 correspondence between 
subspaces of C d and faces of the state space K d . If if is a subspace of 
C d and P is the projection onto H, then the associated face of is 

(31) F P = {A e K d | iuiA C imP} = {AeK d \ imi C 

This correspondence of subspaces of C d and faces of K d follows from, 
e.g., 1 1 . eqn. (3.14)], which says that 

F P = {A e K d | A = PAP}. 
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By (J7j) this is equivalent to ( 13TT) . (Equation (3.14) of [lj is stated in 
terms of positive linear functionals p on associated with the density 
matrices A in via p(X) = tr(AX), but it translates easily to (13T1) 
above.) 

From this it follows that faces of the state space of M m <g> M n are 
themselves "mini state-spaces", i.e., are affinely isomorphic to some K p 
for p < mn. The extreme points of K are precisely the pure states P x , 
where P x denotes the projection onto the span of the unit vector x. 

We recall for use below that the separable state space S is compact, 
as is any face (since faces of closed finite dimensional convex sets are 
always closed.) The extreme points of S are precisely the pure product 
states P X (gty 

We now prove that certain faces of the separable state space are 
themselves "mini separable state spaces", i.e., are affinely isomorphic 
to the separable state space S pq of M p ® M q for some p < m,q < n. 

Notation. If V, W are subspaces of C m , C n respectively with dim V = 
p, dim W = q, then Sep(V" ® W) denotes the separable states in M m <g> 
M n that live on V ®W (i.e., whose image is contained in V®W). Note 
Sep(V r <S> W) is affinely isomorphic to the separable state space S pg . 

We will make frequent use of the following implication for subspaces 
V GC m ,W GC n : 

for x G C m , y G C n , x ® y E V <g> W =^ x £ V and y E W, 

which follows immediately by expanding bases of V and W to bases of 
C m and C n and expressing x and y in terms of these bases. (Alterna- 
tively, cf. [7J eqn. (1.7)]. 

Lemma 9. Let A e M m , B e M n be density matrices. Then 
f&ce s (A ® B) = Sep(imA <g> imB). 

Proof. Note that both sides are compact convex sets, and hence are 
the convex hull of their extreme points. The extreme points of both 
sides will be pure product states, so we can restrict consideration to 
such states. 

Suppose P x ® y € faces(A ® B). This is contained in face^y! ® B), 
which consists of the density matrices whose images are contained in 
im(v4 (g) B) = im A ® imB. Thus x G im.A and ?/ e imB, so P x ^ y G 
Sep(imv4 (g) im5). Thus we shown 

face 5 (A <g> B) C Sep(im A <g> imB). 

For the opposite inclusion, suppose P x ® y is any extreme point of 
Sep(im A g) vmB). Then x ® y G imAtgiimi? implies that x E imA 
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and y G im.B. Hence P x is in face^A) and P y G face^-(P), so there 
exists a scalar A > such that XP X < A and \P y < B. Then 

A® B — [{A - \P X ) + XP X ] ® [(B - XP y ) + XP y ] 

Expanding the right sides gives four separable (unnormalized) states, 
and hence 

P x ®P y = P xm G faces (A ® B). 

Thus 

Sep(im A <g> imB) C face 5 (A <g> B), 
which completes the proof of the lemma. □ 

Lemma 10. Let Vi, . . . ,V q be subspaces of C m and W\, W 2 , ■ ■ ■ , W q 

independent subspaces ofC n . IfO ^ x®y G C m ©C™ ; let J = {7 | x G 
Vy}. Then 



(32) i8|/G^y 7 ®^ 

7=1 

z/f J is nonempty and y G J2-ye.J 
Proof. Assume (152"]) holds. Then 



7 



x®!/e ^ V" 7 <g> W 7 c I v 7 J ® ( ^2 w i 

7=1 V 7 / V 7 

so x G ^7 ^7 an d y € VKy. Thus without loss of generality we may 
assume £ 7 V 7 = C m and ^ W 7 = C n . 

Let Pi, . . . , P q be the (non-self-adjoint) projection maps correspond- 
ing to the linear direct sum decomposition C n = W\ © • • ■ © W q . Then 
for 1 < $ < q 

(33) x® Ppy = {I ® P p ){x®y) eV p ®Wp. 

If we choose j3 so that P^y 7^ 0, then x G Vg, so J is not empty. Then 
for 7 ^ J, we have x ^ V 1 , so by ( l33l) . P 7 y = 0. It follows that 

Conversely, suppose J is nonempty and y G Xl 7 ej ^7; sa y 2/ = 
E 7e jZ/7- Then 

x©y = V'x©y 7 GV'V 7 © W 7 . 



/ , f 'o' yy ^- / j " 7 
7GJ 7=1 



□ 
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We say the convex hull of a collection of convex sets {C a } is a direct 
convex sum if each point x in the convex hull has a unique convex 
decomposition x = X a x a with x a G C a . In the theorem below, 
co denotes the direct convex sum. 

Theorem 11. Let T = A 7 <g> £? 7 be a density matrix in M m <g> M n . 
Assume that A\,...,A P are density matrices with pairwise disjoint 
ranges, and that Bi, . . . , B p are positive matrices with independent im- 
ages. Then the face of the separable state space S generated by T is the 
direct convex sum 

p 

(34) face 5 T = co (J) face s (Ay g> B 7 ). 

7=1 

Proof. We first show that the convex hull on the right side of (1511) 
is a direct convex sum. First note that by the assumption that the 
images of the _B 7 are independent, it follows that the subspaces im Ay® 
imB 7 are independent. (Indeed, combining product bases of imAi (g> 
im B\ , . . . , im A p <g> im B p gives a basis of Yl-y i m ^7 ® im £? 7 , from which 
the independence claim follows.) 

Now suppose C 7 , D n G faces (Ay ® £? 7 ) for 1 < 7 < p, and that 

(35) ^L7 7 = ]TD 7 . 

7 7 

Then for any (GC m ® C n 

7 7 

Since 

faces (Ay <8> 5 7 ) C face^Ay <g> 5 7 ) = {E G if | imi? C im(Ay (g) J B 7 )}, 

then for each 7, C 7 £ and are in im(Ay eg) i5 7 ) = imA 7 ® imi? 7 . 
Hence by independence of the subspaces im A 1 ®im.B 1 , we must have 
C 7 £ = _D 7 £ for each 7 and each vector £. Therefore C 7 = -D 7 for all 7, 
showing that the convex hull is indeed a direct convex sum. 

Next we prove the equality in f[34")) . Suppose P x ® y is in the left side. 
Since the face that the state P x ®y generates in S is contained in the 
face this state generates in K, then x <8> y is contained in the image of 
^ 7 A 7 Cg)£> 7 , which is im A 7 ®im_B 7 (cf. Q). Since by assumption 
A±, . . . , A p are disjoint, the set J in Lemma [10] is a singleton set, so 
there is some such that x G imA^ and y G imB^. Then by Lemma 
El P x ®y G lajces{Ap®Bp), which shows the left side of (1341) is contained 
in the right. 
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The extreme points of the right side are each contained in some 
faces (A^ <g> Bp), and since Ap <g> Bp is one of the summands on the left, 
then 

iace s (Ap <g> Bp) C face 5 (^ A 7 ® B 7 ), 

7 

which completes the proof of fl34|) . 

□ 

Lemma 12. If x is a unit vector in C m and B is a density matrix in 
M n , then 

(36) face 5 (Px ® B) = face^ mn (P, <g> B) = P, <g> face Kn B. 

Proof. By Lemma [HI 

faces(Br <8> B) = Sep(imP x . £g> imB) = Sep(Cx £g> imB). 

Since every vector in Cx ® imB is a product vector, every density 
matrix whose image is contained in Cx ® imB is separable, as can be 
seen from its spectral decomposition. Thus by (I3T1) 

Sep(imP 1 .®imB) = {E e K mn | imB C im(P x ®B)} = face^ m „ {P X ®B), 

so the first equality of (136]) follows. 

Now we prove the second equality of (1361) . If T = P x <g> A with 
A G facex n B, then 

imT = im(P x <g>]m.A) C im(P x . <g> B) 

implies that T G facex m „(Pr <8> B), so we've shown 

P x . <g> face*,, B C face^jP, ® B). 

To prove the reverse inclusion, let T G facex mre (P c <g> B). Then imT C 
Cx <g> imB. We will prove there exists A G face^ n B such that T = 
P x <g> A, which will complete the proof of the lemma. 

Since imT C Cx ® imB, for each y G C n there exists a unique w G 
imB such that T(x®y) — x®w. Define A G M n by x®^4y = T(x®y) 
for y G C n , and observe that imA d imB. 

For z G C m , if z = x or z _L x we have T(z®y) = (P x ® A)(z®y). It 
follows that T = P x ® A. Since T is a density matrix, it follows that A 
also is a density matrix. Since imAc imB, then A G face# n B. Thus 
T eP x ® face Kn B. 

□ 

In Theorem [TTJ the faces of the separable state space are expressed in 
terms of other (smaller) separable state spaces. In some circumstances, 
these are actually state spaces of the full matrix algebras, as we now 
show. (This generalizes [21 Thm. 4].) 
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Theorem 13. Let T = ^7=1 ^7 ® B 1 be a density matrix in M m £g> 
M n . Assume that A±, . . . ,A P are rank one density matrices, and that 
Bi, . . . , B p are positive matrices with independent images. Then there 
are unit vectors x%, . . . x q , with P Xl , . . . , P x distinct, and independent 
density matrices C\, . . . , C q in M m , such that T admits the convex de- 
composition 

q 

(37) T = Y. X - P ^ ® C - 

This decomposition is unique, and the face of S generated by each P Xv £g> 
C v is also a face of K mn , so that 

q q 

(38) face 5 T = co face Xmn {P Xv <8> C v ) = co ®(P Xv ® face^ n C v ). 

v=\ u=l 

Proof. By assumption, each v4 7 is a postive scalar multiple of a projec- 
tion P x , where unit vector in C m . Absorbing this scalar into 
By, we write the given decomposition in the form 

v 

T = Y,P X ~,®B T 
7=1 

Now we collect together terms where the first factors P x coincide. In 
precise terms, we define an equivalence relation on the indices {1, . . . , p] 
by 7 ~ k if Cx 7 = Cx K , or equivalently if P Xl = P Xk . Let J be the set 
of equivalence classes, and for each equivalence class v G J choose a 
representative 7 G 1/ and define Then 

where C u = J2-y& B*t- Define q — \J\ — q; numbering the members of J 
in sequence gives a decomposition of the form specified in the theorem. 

Since the images of the P Xv are disjoint, the final statement of the 
theorem follows from Theorem [11] and Lemma [12j 

□ 

Decompositions into pure product states 

If T is B- independent, Theorem [6] provides a canonical way to decom- 
pose T. Then with the notation of Theorem [6], we can decompose each 
Ay and By further via the spectral theorem into linear combinations of 
rank one projections, and this gives a representation of T as a convex 
combination of pure product states. The next result describes when 
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this decomposition into pure product states is unique, generalizing the 
uniqueness result in [21 Corollary 5]. 

Theorem 14. If T G M m <8> M n is a B -independent density matrix, 
then there is a unique decomposition of T as a convex combination of 
pure product states iff in the canonical decomposition (|22|) of Theorem 
0, each Aj and each B 1 has rank one. Thus the decomposition of T 
into pure product states is unique iff T can be written as a convex 
combination 

p 

7=1 

with unit vectors y%, . . . ,y p that are linearly independent, and unit vec- 
tors x\, . . . , x p such that P Xl , . . . , P Xp are distinct. 

Proof. Suppose that T is 5-independent and admits a unique decom- 
position as a convex combination of pure product states. Let T = 
^7 ® -^7 De the canonical decomposition of T given in Theorem [6j 
If any Ap does not have rank one, then there are infinitely many ways 
to write Ap as a convex combination of pure states, which when com- 
bined with any decomposition into rank one projections for the other 
v4 7 and each 5 7 gives infinitely many decompositions of T into pure 
product states. The same argument applies if any £> 7 does not have 
rank one. Hence if T admits a unique convex decomposition into pure 
product states, each A 7 and B 1 must have rank one. 

Conversely, assume that T can be written as a convex combination 

T = J2 A ^ ® p y-< 

7 

with {P x -y} distinct and with {?/ 7 } independent. This decomposition 
satisfies the hypotheses of Theorem [T3| and thus faces (T) will be the 
direct convex sum of the singleton faces {P Xl ®P Vl }- Now suppose that 
we are given any other convex decomposition into pure product states 

t = y,up Zv ^p Wv , 

where we are not making any assumption about independence of {P Wv } 
or distinctness of {P Zv }- Then each P Zv ® P Wv is in faces T and is an 
extreme point of the separable state space S. By the definition of a 
direct convex sum, we conclude that each P Zv ® P Wv must coincide with 
some P Xi ®P y „- Thus the convex decomposition of T into pure product 
states is unique. □ 
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Remark. One might suspect that for the uniqueness conclusion in The- 
orem [TH it would suffice for the joint eigenspaces of the (T)^ to be 
one dimensional, but this is not correct, as can be seen by considering 
A <g> P y where rank A > 1 . 

The marginal rank condition 

In this section we specialize previous results to an important class of 
separable states. 

Definition. A density matrix T e M m eg) M n satisfies the marginal 
rank condition if rankT = maxfrankT^, rankTg), which reduces to 
rankT = rankTg if m < n, which we will assume in the sequel. 

We will see that such matrices, if separable, are B- independent. The 
following lemma for states on M m <g> M n appeared for m = 2 in [13] . 
and for general m,n in jTUl Lemma 6, and proof of Theorem 1]. An 
alternate shorter proof for general m, n can be found in [TBI Lemma 
13]. It shows that separable density matrices satisfying the marginal 
rank condition are the same as those that admit a representation (pQ) 
with each A 1 and _B 7 of rank one, and with Bi, . . . , B p independent. 

Lemma 15. Let T be separable. Then T admits a decomposition T = 
Y7i=i KPxi® yi with y 1: ...,y p independent iff rankT = rankT B . 

We now show that Theorem Ogives a practical way to check whether 
a particular matrix satisfying the marginal rank condition is separable. 
Theorem [6] then also provides a way to find an explicit representation 
of T as a convex combination of tensor products of positive matrices. 
(For testing separabiltiy, the PPT test also suffices, cf. [TO].) 

Theorem 16. Let T e M m £g> M n with rankT = rankTe. Define T 
as in (j3J). Then T is separable iff the matrices (T)ij are normal and 
commute. 

Proof. Assume T has marginal rank. If T is separable then by Lemma 
[TBI T is B- independent, and hence by Theorem El the matrices (T)^- 
are normal and commute. Conversely, if these matrices are normal and 
commute, then by Theorem[6]T is B- independent, and hence separable. 

□ 

Corollary 17. If T is separable of marginal rank, then T admits a 
unique decomposition 

p 

T = Y J P*- i ®B 1 , 

7=1 
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with Bx, . . . , Bp positive matrices with independent images, and x\, . . . , x. 
unit vectors with P X1 , . . . , P Xp distinct. 

Proof. By Theorem El there is a unique decomposition 

T = A i ® B ~i 

7 

where each A 1 is positive with trace 1, A±,...,A P are distinct, and 
B...,B P are positive and have independent images. By the marginal 
rank condition, 

rank T = rank A^ rank S 7 = rank Tb = rank B^ 

7 7 

which implies that rank A,, = 1 for all 7. Thus there are unit vec- 
tors x 7 such that A^ = P x . Distinctness of A\,...,A P implies that 
P Xl ,.. .,P Xp distinct. □ 

Remark. Uniqueness of the decomposition in Corollary [T7] was first 
proved in [21 Section IV]. Corollary [17] provides an alternate proof, and 
Theorem [HI provides an explicit way to find that decomposition. 

Summary 

We have defined a class of separable states (5-independent states) 
which generalizes the separable states whose rank equals their marginal 
rank. We have described an intrinsic way to check membership in this 
class, and for such states we have given a procedure leading to a unique 
canonical separable decomposition into product states. 
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